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G\ Abstract: The moment of inertia of a trapped superfluid gas of atomic Fermions 

^ ( 6 Li) is calculated as a function of two system parameters: temperature and 

C deformation of the trap. For moderate deformations at zero temperature the moment 

>— > of inertia takes on the irrotational flow value. Only for T very close to T c rigid rotation 

(N| is attained. For very strong trap deformations the moment of inertia approaches its 

^ rigid body value even in the superfluid state. It is proposed that future measurements 

^ of the rotational energy will unambiguously reveal whether the system is in a 

q superfluid state or not. 

On 

S 1. Introduction 

a 

q The advent in 1995 of Bose-Einstein-Condensation (BEC) of atomic Bosons in 

magnetic traps certainly represents a milestone in the study of bosonic many 
body quantum systems. This is so because a systematic study of these systems, 
starting with the free particle case, as a function of increasing density, particle 
number, and other system parameters seems possible and has largely already 
been cut into and is progressing with a rapid pace [1,2]. On the other hand the 
recent experimental achievement of trapping 6 Li atoms and other fermionic alkali 
atoms [3] also spurs the hope that for the fermionic many body problem, as much 
progress will be made in the near future as for the bosonic systems. For instance 
the fact that trapped spin-polarized 6 Li atoms which feel a strong attraction in the 
triplet channel (scattering length: a = -2063 a with a the Bohr radius) may 
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undergo a phase transition to the superfluid state has recently provoked a number 
of theoretical investigations [4]. One major question which is under debate is how 
to detect the superfluidity of such a fermionic system, since in contrast to a 
bosonic system the density of a Fermionic system is hardly affected by the 
transition to the superfluid state [5]. Several proposals such as the study of the 
decay rate of the gas or of the scattering of atoms off the gas have been 
advanced [4]. Though such investigations may give precious indications of a 
possible superfluid phase, we think that in analogy with nuclear physics, a 
measurement of the moment of inertia certainly would establish an unambiguous 
signature of superfluidity. To measure the spin and the rotational energy of 
trapped atoms definitely is a great challenge for the future. However, in nuclear 
physics, where y-spectroscopy is extremely well developed, the strong reduction 
of the moment of inertia with respect to its rigid body value has been considered 
as a firm indicator of nucleon superfluidity immediately after the discovery of 
nuclear rotational states almost half a century back [5]. Therefore awaiting future 
experimental achievement also for trapped fermionic atoms, it is our intention in 
this work to give some theoretical estimates of the moment of inertia as a function 
of some system parameters such as deformation of the traps or temperature of 
the gas. In this study we can largely profit from the experience nuclear physicists 
have accumulated over the last decades in describing such phenomena. The 
expectation is indeed that there will be a great analogy between the physics of 
confined atomic Fermions and what one calls in nuclear physics the liquid drop 
part of the nucleus. As astonishing as it may seem assemblies of fermions 
containing no more than -200 particles (nucleons) already exhibit an underlying 
macroscopic structure well known from the Bethe-Weizsaecker formula for 
nuclear masses [5]. In superfluid rotating nuclei as early as 1959 Migdal proposed 
a statistical description of the nuclear moment of inertia [6] which grasped the 
essential physics of a self contained rotating superfluid Fermi liquid drop and 
which serves as a reference even today. 

In the present work we will cast Migdal's approach into the more systematic 
language of the Thomas Fermi theory which together with its extensions is 
applied extensively since decades to normal fluid but also to superfluid nuclei 
[5,7,8]. It is fortunate that we can profit from this experience for the description of 
trapped fermions, since their number of order 10 5 , together with the smoothness 
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of the potential, certainly turns a statistical description into a very precise tool. On 
the other hand it may not be excluded that in the future the study of much smaller 
systems of trapped atomic Fermions with numbers ~10 2 may be studied probably 
revealing many analogies with real nuclei such as shell structure etc. The 
investigation of the transition from microscopic to macroscopic as the number of 
particles is increased continuously may then become a very interesting field also 
in the case of atomic Fermions. In detail our paper is organized as follow. In 
sections 2 and 3 we review the Thomas Fermi approach to inhomogeneous 
superfluid Fermi systems. In section 4 first the so-called Inglis part of the moment 
of inertia of a rotating superfluid and confined gas of atomic Fermions is 
presented. Second the influence of the reaction of the pair field on the moment of 
inertia is calculated. It is shown that this leads to the irrotational flow value in the 
limit of strong pairing. In section 5 the current distributions in the superfluid and 
normal fluid regimes are contrasted. In section 6 the numerical results are 
presented in detail. Finally in section 7 we discuss our results and present some 
conclusions. 

2. Thomas Fermi approach to fermionic atoms in deformed traps 

As in the boson case the Thomas Fermi (TF) approach [5] to trapped atomic 
gases is rendered extremely simple by two facts: the smoothness of the traps 
(harmonic oscillator) and the large interparticle distance which makes a 
pseudopotential approximation to atomic interactions valid. Let us therefore write 
down the TF equation for a double spin polarized system of trapped ( 6 Li) atoms 
in the normal fluid state. For convenience we first consider the system at zero 
temperature T discussing the T * case later on. In TF approximation the 
distribution function is given by, 

f(R,p) = 6(ju-H cl ) (2.1) 
with, 
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H cl =^ + V ex (R)-gp(R) 



(2.2) 



Here ju is the chemical potential, and V ex (R) the trap potential supposed to be of 
harmonic form. The density p(R) is obtained from the self-consistency equation, 



the local Fermi momentum. The coupling constant g is related to the scattering 
length in the same way as in the Bose gases [1 ,2 ] via: 



The TF equation (2.3) leads to a cubic equation for the self-consistent density, 
which can be solved without problem as a function of the external potential. In 
this paper our main interest will be the study of the moment of inertia of a rotating 
condensate. Since the study is very much simplified assuming that the self 
consistent potential is again a harmonic oscillator and since the effect of the 
attractive interaction between the atoms essentially results in a narrowing of the 
self consistent potential with respect to the external one we will use instead of the 
exact TF solution for the density the following trial ansatz for the local Fermi 
momentum: 




(2.3) 



with, 




(2.4) 
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kf l (R)- 



\h l 2 



(2.6) 



where 6^ and <a z are the variational parameters. The chemical potential is 
determined from the particle number condition, 



N = \d'rp ,rial {r) 



(2.7) 



and the kinetic energy density is given by, 



7(R)= ^- /(R ' P) = I^^ (R)]5 



(2.8) 



We then can analytically calculate the total energy, 



E(a) ± ,0) z )=^d 3 R 



r(R)+K,(R)p(R)-fp 2 (R) 



(2.9) 



as a function of 6^ and a> z . Minimizing this expression with respect to &i ± and 
&) for a given external deformed harmonic oscillator potential, 



(2.10) 



leads to the variational solution. For the spherical case q) ± =q) z =q), this is 
shown in Fig. 1. We see that this approximation to the TF equation is quite 
reasonable. For co Q = 6.9 s 1 , the value of the variational frequency is co = 7.69 s 1 
that is co>co Q , implying a compression of the density. Increasing a> by 6% from 

its variational value allows an almost perfect reproduction of the full TF solution. 
We will adopt this latter value in all our forthcoming calculations. In order to obtain 
the deformed case we simply squeeze the potential using a deformation 
depending on one parameter, 
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CO, 



(2.11) 



Frequencies co L and co z are then defined as functions of the deformation 
parameter by, 



co L = coS~ V3 ' 
a = coS^ 



(2.12) 



in order to keep constant for any value of 8 the volume, 



(coJco z = 0) 3 



(2.13) 



From now on we therefore will use for the non superfluid Wigner function at zero 
temperature the expression, 



/(R,p) = e{» - f - ^ (col fc +Rl )+ co\R\ f 
2m 2 



(2.14) 



with co z ,co ± from (2.12) and ju determined from the particle number condition. 



3. The superfluid case 



Since trapped spin polarized 6 Li atoms feel a strong attractive interaction in 
the triplet channel the system very likely will undergo a transition to the superfluid 
state at some critical temperature T c as was discussed in detail in ref. [4]. As we 
have pointed out in the introduction the superfluid state will unambiguously reveal 
itself in its value of the moment of inertia. At the moment the measurement of 
angular momenta of trapped Bose or Fermi gases has not been achieved and 
represents a great future challenge to the experimenters. In order to establish 
how the two essential system parameters which are the value of the gap, i.e. the 
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temperature and the deformation of the external trap influence the value of the 
moment of inertia, we will now proceed to its evaluation in the superfluid state. 

Since we are dealing with an inhomogeneous system, even in the non 
rotating case the gap is actually a non local quantity A(r,r') or in Wigner space 
A(R,p). It has been shown in [8,5] that to lowest order in //the gap equation is 
given by, 



. x r d 3 k , ,xA(R,k) J #(R,k)^ 
A(R,p)=J - ^ v(p-k) V a anh 



2£(R,k) 



IT 

\ J 



(3-1) 



where £(R,p) is the quasiparticle energy, 



E(R,p)= 



2m* (R) 



+ A 2 (R,k) 



(3-2) 



with p F (R)=hk F (R) the local Fermi momentum (2.4). Since the effective mass 
m* is so far unknown for the trapped gases of atomic Fermions we will take 
m*=m. Furthermore, for the time being, as in [4], we will eliminate the 
interatomic potential v in (3.1), expressing it by the scattering length (2.5). We 
then obtain [4], 



tanh 



£(R,k) 



IT 



2£(R,k) 2(e k -£ F (R)) 



A(R,k) 



(3-3) 



where P stands for principal value, e k = 



ft k , fi k p . . 
and e F = — - . At zero 

2m 2m 

temperature, as described in [9], (3.3) can be solved analytically in the limit 

A(R,p F (R)) 



e F (R) 



-> . The result is given by, 



7 



A,(R) = A(R,fc F (R))=8e- 2 ^(R)e 2 ^ R) l fl l (3.4) 



A posteriori one can verify that — «1 for all values of R and therefore (3.4) is 

e F 

an excellent approximation to (3.3). This also has been verified in [4]. For 2.865 x 
10 5 6 Li atoms, the case considered in [ 4], the gap is shown for a spherical trap 
as a function of the radius in Fig. 2. 

For the determination of the critical temperature T c and, later on, for the moment 
of inertia we will need the value A of the gap averaged over the states |n)at the 

Fermi level, 

A(£ F hA = Tr{Ap(£ F )) (3.5) 



with, 



P(e F )=-}-. T\n){n\S{e F -e n )=^—S{e F -H) (3.6) 



where |n)and e n are the shell model states and energies and, 

d{s F -£ n )=Trd{s F -H) (3.7) 

n 

is the level density at the Fermi energy. 

It has been shown in [10] that again the TF approximation leads to an 
excellent average value, 

^7k-/Sf AM '- Hj (3 ' 8) 
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In the spherical case with A f (r) from (3.4) all integrals but the radial one can be 
performed analytically, the latter being done numerically. For the case shown in 
Fig. 2 one obtains, 



A = 16.4 nK (3.9) 

One now can discuss the question whether it is the value (3.9) which determines 
the critical temperature T c or whether it is A F (R = 0)~ 65 nK as determined from 
Fig. 2. In the strict quantum mechanical sense the BCS equations should be 
solved in the shell model or better HF basis and then T c is a global parameter 
which must be determined from the solution of the quantum mechanical gap 
equation. Since we believe that the value in (3.9) comes rather close to the 
quantum mechanical value of the gap at the Fermi energy we think that (3.9) 
should determine T c which we can obtain from the usual BCS weak coupling 
relation [11] A = 1.767; to be 

7>10nK. (3.10) 

On the other hand the temperature dependent BCS approach is certainly not an 
exact theory for finite systems and it may be that e.g. fluctuations and 
correlations break off the quantum coherence of the global single particle wave 
function and more local cluster structures are formed. In such cases it could be 
that superfluidity persists locally to higher temperature than the one deduced 
from (3.9). A precise answer to this question is certainly quite difficult and in view 
of this uncertainty we will later for the moment of inertia take the critical 
temperature as a parameter which can vary within certain limits determined by 
the gap parameter (3.9) or A F (o) respectively. Independent of the question of the 
precise value of T c we will have to know the detailed T-dependence of the gap 
A(r) which however, in BCS theory, given A and T c , is determined by the 
a(t) T 

universal function — j-^\r\ terms of — . This function is determined from the 

A(0) T e 

solution of the equation [11], 
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^A(T))_ fA(T) 



-In 



v A(0) y 



T 

V 1 J 



with, 



(3-11) 



A{u)= \dy-F=^ 



o ^jy 2 +u^ 



1-tanh 



For completeness it is again shown in fig. 3. This T-dependence of the gap we 
will later use for the evaluation of the moment of inertia. 

4. Moment of inertia 

The moment of inertia of a rotating nucleus has fully been formulated in linear 
response theory (i.e. RPA) by Thouless and Valatin [12]. The corresponding 
expression is therefore called, in the nuclear physics literature [5], the 
Thouless-Valatin moment of inertia. It consists of two parts, the so-called Inglis 
term, which describes the free gas response, and the part, which accounts for 
the reaction of the mean field and pair potential to the rotation. In the 
superfluid case the Inglis part has been generalized by Belyaev [13] and the 
linear reaction of the gap parameter onto the value of the moment of inertia 
was first evaluated, together with the Inglis term, by Migdal [6]. The reaction of 
the HF field on the rotation is a minor effect and we will neglect it in this work. 
We therefore will write the moment of inertia as a sum of the Inglis-Belyaev 
term ®,_ B and the Migdal term ® M . In total, 



= 7 _ B +0 



u 



(4.1) 



In order to derive an expression for in linear response theory we will use 
the Gorkov approach described in detail in many text books (in what follows 
we will use the notation of [14]). Since in addition the derivation of the linear 
response for is given in the original article of Migdal [6] and rerepresented 
in a more elaborate version in [7], we will be very short here and only give 
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more details where in our opinion the presentations in [6,7 ] may not be 
entirely explicit. Let us start writing down the Gorkov equations in matrix 
notation, 



_ — -H + ju 

dT 



G = l-AF' 



(4.2) 



dT 



F + =A*G 



with 



H = H Q -nL x =H + H l 



(4.3) 



where H is the shell model Hamiltonian (2.2 ) and 

L x= r yP Z - r zPy 

the angular momentum operator corresponding to a rotation with angular 
frequency CI around the x-axis. In (4.2) G and F are the normal and anomal 
Matsubara Green's functions (see Chapt. 51 of [ 14]), 



G nn ,=-(T T a n (T)a/(Tj 



(4.4) 



F\ n >=-(T T a;(T)a/(4 



Linearising (4.2) with respect to H x , that is G = G + G { , F + =F +F X and 
A = A + A, (as mentioned we will neglect the influence of the rotational field on 
H ) one obtains for (4.2), 



(4.5) 



with, 



n 



G U _ B - GqH x Gq + F H 1 F Q 

G IM = _G (A^0 + _ ^0 +A l G 



(4.6) 



and, 



F, + = D H;F Q + + F G + Hfi, - F + A^o + + D A/G, 



o 



(4.7) 



where, 



D = 



O) n 2 +H 2 +A Q 2 



(O n 2 +H Q 2 +A Q 2 



0) n 2 +H Q 2 +A Q 2 



and a> n are the Matsubara frequencies [14]. 

In (4.5,4.6) we have split the first order Green's function in an obvious notation 
into the Inglis-Belyaev and Migdal contributions. For the latter one needs the 
linear reaction of the pair field to the rotation. We later will see how this can be 
determined from (4.7). First let us, however, evaluate the l-B part of the 
moment of inertia. 

4.1 The Inglis-Belyaev part of the moment of inertia 

The l-B part of the moment of inertia can be evaluated without the 
knowledge of A l i.e. without the use of (4.7). The density response 

corresponding to G U _ B of (4.5) is evaluated from the limit z'^>z + or from 
summing over the Matsubara frequencies in the upper half plane (see Ch. 7 of 
[14]) . One obtains the well known result [5, 6, 7, 14], 




(4.8) 
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with, 



F=F + {i-f-r)+F_{f-r) 



(4.9) 



where, 



F ± =F ± (e n ,e n ,) = 



E n E n , + 4 n 4 n ,-A(e n )A(e n ) 
2E H E H .fa±E H .) 



(4.10) 



f = /(£„)= 



1 



f'=f(e n ') 



(4.11) 




and, 



E n =^ n 2 +A 2 kJ 



(4.12) 



are the quasiparticle energies with as ingredient s n , the shell model energies. 
The gap parameters A n have been replaced in (4.10), in analogy to (3.5), to 
statistical accuracy by A(e n ), the ones averaged over the energy shell. The 
moment of inertia is given by , 



Since we are interested at temperatures T<T c , which are very low with 

respect to the Fermi energy, we checked that one can to very good accuracy 
neglect in (4.9) the thermal factors (4.11). The only important temperature 
dependence of the moment of inertia therefore exists via the T-dependence of 
the gap. We thus will henceforth treat all formulas as in the T=0 limit keeping, 
however, the T-dependence of the gap. With this in mind we can write for the 
moment of inertia, 



®i-B=Tr(hPu-n)- 



(4.13) 
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®j_ B = ^ JJ d(odaf8{fo - e n ]S(a)-e n ,)\(n \L X | h')| 2 F + (co, 0)') 



(4.14) 



In this formula the important quantity to calculate to statistical accuracy is, 

L x 2 (n, n') = \(n \L x \nf = Tr[(L J( | n)(n%\n){n |)] 



= \^^Xi\n){n%\n){n\\ 



(4.15) 



where O w = <9(R,p) means the Wigner transform of the operator O [5]. To this 
purpose we again replace the density matrices \n)(n\ and |n')(n'| by their 
average on the energy shell (3.6), 

\n){n\^p(s n ) 

We therefore obtain, 



7 _ B = \\dcoda> '}^/ [(L x \ (L x (co, co% ]f + (<», ©') 



{mi y 



with, 



(L>,<y')) w = \s(p'-H )t x S((Q-H Q )l 



(4.16) 



(4.17) 



Introducing into (4.17) the Fourier representations of the two 6 -functions and 
transforming to center of mass and relative coordinates one obtains, 



with, 



(2;th f 









-iH„T j 




-iH„T 

e 


v 2 ; 





w 



r cq + cq' 

E = ; e = co-co 



(4.18) 



(4.19) 
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and, 

O(t)=e iH °'O(0)e- iH °' 



(4.20) 



To lowest order in ft we replace the triple operator product in (4.18) by the 
product of their Wigner transforms [5], 



lim 

fi->0 

















-m„T 
e 


— o ~ i2H 0cl T ' jcl 

w~ e u x 






v 2 ; 






V 2 , 



and therefore, 

(L, (co, w% = L x (E, e,R,p)=±5(E- H 0cl )\-^^L 



1 T f 

cl 



with 

( 

j cl 
^x 

V 



R, 



) 



-V r 

\ 2 ) \ 2 J 



-R 



(4.21) 



(4.22) 



(4.23) 



At this point the choice of our approximate self consistent potential of 
harmonic oscillator form (see 2.14) turns out to be very helpful, since the 
classical trajectories in (4.23) can be given analytically. 



R, 



Pi 



v 2 / 



v 2 y 



R; COS 



V 



'2 



Pi 

H ' — COS 

mCO, 



hco. 



'2 



(4.24) 



= p,cos 



flCO: — 

2 



f 



+ mO) ; R ; cos 



hco, 



V 



'2 



with i = x,y,z and a) x =d) y =d) ± 

In the phase space integral of (4.15), for reasons of symmetry, only the 
diagonal terms of L c x .L c x i^) survive and therefore we obtain, 



cos 



flCO, 



f 



+ 









sin tiCQ 1 — 


sin 


hco, - 


I I 2 J 




v ~ 2 ) 



(4.25) 



where, 
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.TF 



1 

6^ 



2m 



-i 3/ 



(E-V) 



(4.26) 



is the density in TF approximation. 

The product of cosine and sine in (4.25) can be expressed in terms of the 
cosine of the sum and difference of the arguments and then the t -integral in 
(4.23) can be performed. This leads to S functions what allows doing also the 
e -integral. Furthermore, as shown by Migdal [6], 



F + (E,eh 



1-G 



v 2A y _. 



8{E-n) 



where (see eq. (3.5)), 



A = A(e F ) 



(4.27) 



and, 



G(x) = 



arsinh(jc) 
x^l + x 2 



(4.28) 



Finally one obtains for the l-B part of the moment of inertia the following 
analytical expression [ 6,7], 



= 



rigid 



G + Q)_ 2 + G_Q) + 2 

2 2 

CO 



(4.29) 



where, 

co ± =co ± ±co z ; G ± = G 

and, 



' hco ± 



_ co z +co ± 
rigid co z 3 co ± 4 



12ft 3 , 

v J 



(4.30) 



(4.31) 



is the moment of inertia of rigid rotation. From (4.29) we see that 
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limft = (?) ■ lim (?) =0 

A->0 A->°<= 



The latter result is clearly unphysical and we will see how the account of the 
reaction of the pair field on the rotation will reestablish the physical situation. 

4.2 The Migdal term 

The density response corresponding to the Migdal term is obtained from 
(4-6), 

In (4.32) we need to know A l which we can gain from (4.7) in the following 
way; in the limit i ->? + we obtain from F^the anomal density , 

2£ nJ B„,(£„ +J B„,) 



In analogy with the non-rotating case where tc = — , we also have, 

22s 



A 1 1 A 

k) .=-a- — +— 



(4.34) 



This relation stems from the fact that the quasiparticle energies contain the 
gap only in the form AA* and therefore there is no further first order correction, 
since in our case the external field is a time odd operator and thus, 

A\=-A^-iCl X (4-35) 

Equating (4.33) and (4.34) yields, 
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i«n'A „' + (a o „ + A 0n . 

? r ^ = (4.36) 

2E E ,(E +E ,) 

At this point we again exploit the fact that expression (4.36) is strongly peaked 
around the Fermi energy surface. Following [6], in analogy with (4.27), we 
have, 



[E n E n ,(E n +E n ,)Y 



1 



2A 



£ +£ , 



With (4.35) we then obtain for (4.36), 



(4.37) 



(n\L x \n , J ( e_ —e_, ^ 



2A 



X n 



£ ~£ , 

n n 



2A 



£ +£ , 

n n 



= 



(4.38) 



where L x stands for the time derivative ofL c . Summing on n and n' and following 

exactly the same line of semi-classical approximations as the ones used for the 
derivation of 7 _ B one arrives at the following relation [7], 



\dTG{x)[ 



d 3 p 



(2^z) 3 [2A 4A 2 _ 



S(ju-H 0cl )=0 



(4.39) 



where G{x) is the Fourier transform of G{x) (4.28). 
For the potential in (2.14) (4.39) is solved by, 

%(R)=aR y R z 

with, 

G^ + G 



a = -2Ameo+eo_ 



0) 2 + G + + w 2 _G_ 



(4.40) 



(4.41) 



Inserting this solution into (4.32) leads for the Migdal part of the moment of inertia to 
[6,7], 



W /W ^ rigid 



(olcol (G + + G_f 
col + °>- °>l G + + M-G- 



(4.42) 



Together with (4.29) the expression for the moment of inertia is now complete. Let us 
again mention that we neglected the temperature dependence besides the one 
contained in A = A(r), since all other T-dependence for T <T C is negligible. The 



18 



moment of inertia can then be calculated as a function of deformation and 
temperature. For example it is immediately verified that for A -> oo (4.42) yields the 
irrotational flow value, 




(4.43) 



and therefore, 



lim © = lim + M ) = 



irrot 



(4.44) 



what is the correct physical result. 



5. Current distribution 

One other quantity, which may be interesting also from the experimental point 
of view, are the current distributions of the superfluid rotating gas. Indeed after a 
sudden switch off of the (rotating) trap the atoms will expand keeping memory of their 
rotational state. So if the velocity distribution of the expanding atoms can be 
measured, one may be able to deduce the rotational motion the atoms have had 
before the trap was taken away. The current distribution, as we will see, depends, as 
the moment of inertia, strongly on the superfluid state of the gas. In order to calculate 
the current distribution we first write down the Wigner function of the density 
response which can easily be red off from the formulas given in Section 4. In obvious 
notations we obtain [7], 



A/ _ B (R,p)=^[Rxpl^-// J 



(5.1) 



-a 



Q) + G_ - CO_G 
2^ 



+ 
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With the usual definition of the current, 
one obtains, 



(5.3) 
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(5.5a) 
(5.5b) 



with of course, j x =0. 

Again we see that in the limit A -> oo the current approaches the correct irrotational 
flow limit, 



j^-2 fe n«v( r/l ) 

a-^ co L + co z 



(5.6) 



whereas in the limit of A -> Owe obtain a rigid body current. As we have seen for © , 
as a function of temperature and deformation, we easily can go from one limit to the 
other. We therefore show in Fig. 4 a,b the current distribution for the two extreme 
cases of irrotational and rigid body flow in the laboratory frame. We see that the flow 
pattern is completely different in the two cases. As a function of temperature one 
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continuously can pass from one flow pattern to the other. The point we want to make 
is that for small deformations, as can be seen from (5.6) there is almost no 
irrotational current for low temperatures and this should be detectable experimentally. 

6. Results 

In Fig. 5 we show © as a function of A(r) (and with Fig. 3 as a function of T) for 

several values of the deformation parameter <5 _1 . We see that, not unexpectedly, for 
T ~0 the system is in the irrotational flow limit. This means that the number of levels 
in a range A(o) is very high and thus our statistical assumptions well verified. Only for 
temperatures very close to the critical temperature where A<10nKthere is then a 
strong turn over to the rigid moment of inertia. We also see that the details depend in 
an appreciable way on the deformation: the greater S the earlier starts the bending 
to the rigid value. Tracing the moment of inertia as a function of deformation for 
different values of the gap yields the curves shown in Fig. 5. We see that, 
independently of the gap value, for sufficiently strong deformation the moment of 
inertia will always end up in its rigid body value. This is easily understandable since in 
the limit of very strong prolate deformation the system becomes essentially one 
dimensional. Turning a fluid contained in a rod shaped trap, whether superfluid or 
not, one always has to turn the whole body of the fluid equivalent to rigid body 
rotation. We therefore see that the moment of inertia of a trapped superfluid gas of 
atomic fermions can span a very wide range of values going from almost zero at low 
temperatures and low deformations of the trap to the rigid body value at temperatures 
close to T c and/or at strong deformations. Once the problem of putting the trapped 

gas into rotation is mastered experimentally and the question of measuring the 
moment of inertia solved one will then be able to demonstrate very clearly in 
changing the system parameters whether or not the system is in a superfluid state. 

7. Discussion and conclusions 

In this work we transcribed work of nuclear physics for the evaluation of the 
moment of inertia of a rotating superfluid nucleus to the situation of a rotating gas of 
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trapped superfluid atomic fermions. The question of how to experimentally detect the 
superfluid state of such a gas has recently been discussed in the literature to quite 
some extent. We argue that, like for nuclei, one unambiguous indication for 
superfluidity can come from the study of the moment of inertia as a function of such 
system parameters as temperature and trap deformation. We have seen that the 
value of can run through large variations as the gap, as a function of T, varies 
from zero to its maximal value. The same is true for the variation of © with the 
deformation of the trap. The direct determination of the moment of inertia may be 
difficult. However the variation at constant deformation of the rotational energy, 




(7.1) 



as a function of the temperature, i.e. as a function of A(r) should be measurable 
scanning the range from & rigid > > & irrot . Indeed taking for example a rotational 

frequency of Q = 1.6s _1 (which should be compared to our harmonic oscillator 

frequency a> ~7s~ l ), we obtain for S = co y^ ) =0.5176 the following variation (see 

Fig. 6) for , 



e = 



(7.2) 



as a function of a(t). 

We see that the rotational energy makes about 30% of the harmonic oscillator energy 
for rigid rotation, dropping by a factor of ~3 for superfluid rotation. Supposing an 
experimental error of -10% in the measurement of the rotational energy obtainable 
from a sudden switch off of the trap, one deduces that the variation as a function of T 
of the rotational energy should be in the range of experimental accuracy. However, 
other ways to detect the rotational dynamics may be imagined. Suppose an almost 
spherical rotating trap is suddenly stopped from rotating. For temperatures T>T C there 
is rigid rotation and the Fermi gas will continue rotating for a while. On the other hand 
for T«T C the trap will rotate but the superfluid gas will not be in a rotational state 
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neither before nor after the stop of the trap. Whether the gas rotates or not may 
eventually be detectable by optical means. 

In above discussion we ignored the possibility of vortex formation. The onset of 
instability versus vortex formation in a finite Fermi system is not a completely easy 
task and we will postpone such an investigation to a future work. However, since the 
rotational frequencies Q. considered in the present work are much smaller than the 
oscillator constant a> {&/ < 0.2) we may hope that our results will not be spoiled by 

the appearance of vortices. An indication can also come from the case of trapped 
Bosons where vortices depending somewhat on the number of atoms, do not appear 
for values a/ < 0.5 (see ref.[2]). 

From the above discussion we see that it may well be in experimental reach to reveal 
an eventual superfluid state of the gas once the technique of putting the trap into 
rotation will be mastered experimentally. 

In short we proposed in this work to measure the dynamics of a rotating trapped gas 
of atomic Fermions as a function of temperature and deformation to detect whether 
the system is in a superfluid state or not. Quite detailed and quantitative calculations 
for the moment of inertia and velocity distributions have been presented. Other 
quantities well studied in the case of rotating superfluid nuclei [5] such as Yrast lines, 
even-odd effects, particle alignment, etc., may also become of interest in this case. 
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Figure captions 



Figure 1 : Density profiles ,for the case of a spherical trap, of the non-interacting 
case (full line), the interacting case once calculated exactly from (2.3) (crosses) 
with V ex given by (2.10) and once using the variational^ determined harmonic 

oscillator potential(open squares). Squeezing the variational co by 6% yields a 
density which lies on top of the exact TF solution. 

Figure 2 : The gap for a spherical trap as a function of the radius 

Figure 3 : Ratio of the energy gap to the gap at T=0K as a fun ction of 
temperature 

Figure 4 a and b : The current distribution for the two extreme cases of 
irrotational and rigid body flow in the laboratory frame (note the change of scale 
between the two figures) 



Figure 5 a and b :The moment of inertia as a function of deformation for different 
values of the gap (a) and as a function of the gap for different values of the 
deformation (b) 

Figure 6 : The variation of e (eq. 7.2), ratio of the rotational energy of the gas 
and the harmonic oscillator energy, as a function of the gap energy A(T) for 

different values of deformation S = — and rotational frequency Q.. 
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